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ABSTRACT 

The nearby brown dwarf LP 944—20 has been monitored for radial velocity 
variability at optical and near- infrared wavelengths using the VLT/UVES and 
the Keck/NIRSPEC spectrographs, respectively. The UVES radial velocity data 
obtained over 14 nights spanning a baseline of 841 days shows significant vari- 
ability with an amplitude of 3.5 kms" 1 . The periodogram analysis of the UVES 
data indicates a possible period between 2.5 hours and 3.7 hours, which is likely 
due to the rotation of the brown dwarf. However, the NIRSPEC data obtained 
over 6 nights shows an rms dispersion of only 0.36 kms -1 and do not follow the 
periodic trend. These results indicate that the variability seen with UVES is 
likely to be due to rotationally modulated inhomogeneous surface features. We 
suggest that future planet searches around very low-mass stars and brown dwarfs 
using radial velocities will be better conducted in the near-infrared than in the 
optical. 

Subject headings: stars: low-mass, brown dwarfs 

1. Introduction 

Brown dwarfs (hereafter BDs) are objects with very low-masses (M < 0.075 M & ) that 
develop degenerate cores before being able to settle on the stellar main-sequence (Kumar 1963) 
Large numbers of BDs have been identified in the field, in open clusters and in star- 
forming regions. These substellar-mass objects could be as numerous as the known stars 
(Bejar et al. 2001; Chabrier 2003). 
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The detection of strong emission lines in young BDs, as well as near-infrared, mid- 
infrared and millimiter emission exceeding those expected for normal photospheres, strongly 
indicate that young BDs have disks where planets may form (Barrado y Navascues & Martin 2003; 
Klein et al. 2003; Luhman et al. 2005). Moreover, mid-infrared spectroscopic observations 
obtained with Spitzer/IRS indicate dust evolution in the discs around young BDs (Apai et al. 2005). 

One of the first obvious targets to search for planets around BDs is LP 944—20 (BRI B0337- 
3535; LEHPM 3451; 2MASS J03393521-3525440; J 2M ass =10.72, ( J - K) 2MASS =1.18) be- 
cause it is nearby (d=5.0 pc, (Tinney 1998)) and young (age~320 Ma, (Ribas 2003)). This 
BD was first spotted over 30 years ago (Luyten & Kowal 1975) as a dim high proper- motion 
red star, and its substellar nature was revealed through the spectroscopic detection of lithium, 
a trademark of BDs (Magazzu et al. 1993; Tinney 1998). LP 944— 20 has a relatively early 
spectral class for a BD (dM9, (Martin et al. 1999)) because of its young age, and has a low- 
level of chromospheric and X-ray coronal activity (Martin & Bouy 2002), although occas- 
sional flares have been reported (Rutledge et al. 2000). Nevertheless, LP944-20 was detected 
as a radio source suggesting the presence of strong magnetic activity (Berger et al. 2001). 

In this paper we present a search for giant planets around LP 944—20. We have moni- 
tored the radial velocity (hereafter RV) of this BD over 20 nights at optical and near-infrared 
wavelengths. We find that the RV data obtained in the optical is variable with an amplitude 
of 3.5 kms -1 and a period between 2.5 and 3.7 hours. However the RV data measured in our 
near-infrared spectra have an rms dispersion of 0.36 kms" 1 and does not follow the periodic 
pattern seen in the optical. We discuss the implications of these results for RV searches for 
planets around brown dwarfs. 

2. Observations and results 

The visible spectroscopic observations of LP 944—20 were obtained with the Uv- Visual 
Echelle Spectrograph (UVES) on the VLT Unit telescope 2 (KUEYEN) in service mode. 
The instrumental configuration was the same as described in (Guenther & Wuchterl 2003). 
This setting covers simultaneously the regions from 667.0 nm through 854.5 nm and from 
864.0 nm to 1040.0 nm. It includes the telluric bands between 686.0 nm and 693.0 nm and 
between 760.0 nm and 770.0 nm. The telluric lines were used as secondary wavelength refer- 
ence because they are know to be stable within 15 ms" 1 (Balthasar et al. 1982; Smith 1982; 
Caccin et al. 1985). Six UVES RV measurements of LP 944—20 have been reported in the lit- 
erature and a possible RV variability of this object has been noted (Guenther & Wuchterl 2003). 
In this paper we report nine additional UVES RV points. All of them were computed in the 
same way as described in detail in previous work (Guenther & Wuchterl 2003). 
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Our infrared spectroscopic observations of LP 944—20 were obtained with NIRSPEC on 
the Keck II telescope. They have already been described in the context of the measurement 
of rotational broadening (Zapatero Osorio et al. 2006). We obtained one more NIRSPEC 
RV datapoint in January 2006 which was not included in the rotational velocity paper. The 
data reduction was identical to that of the previous spectra. Near-infrared measurements 
of the heliocentric RV of LP 944—20 were computed via cross correlation of the NIRSPEC 
data with spectra of the field vB 10 (GJ752B, dM8), which was observed with very similar 
instrumental configuration in 2001 June 15 and 2001 November 02 (see Fig. 1). From our 
NIRSPEC data, we measured the radial velocity of vBlO to be +34.7 ±1 kms" 1 , in very 
good agreement with the various values (+35.4, +35.2, +35.0 kms -1 ) found in the literature 
(Martin et al. 1996; Mohanty & Basri 2003). We adopted +35.0 kms -1 for the correlation 
procedure. In addition, vB 10 is known to be one of the slowest rotators among ultracool 
dwarfs (wsini = 6.5 kms" 1 , (Mohanty & Basri 2003)). 

We calculated the cross-correlation function in multiple echelle orders of the NIRSPEC 
spectra using the task FXCOR within iraf, and measured velocities by fitting a gaussian to 
the cross-correlation peak. We used 4 spectral orders that are free from strong telluric line 
contamination. Two examples of those orders are shown in Fig. 1. These orders include 
features due to Kl, TiO, FeH, and water vapor. We note that given the good signal-to-noise 
ratio of the spectra and the significant similarity in spectral type between LP 944—20 and 
vB 10, the peak of the cross-correlation function is unambiguously identified. The NIRSPEC 
radial velocities of LP 944—20 shown in Table 1 are the mean of the results in the different 
orders. The associated error bars correspond to the standard deviation of the mean. 



3. Discussion 

The UVES RV dataset shows a significant scatter and we performed a periodogram 
analysis of it. The periodogram analysis using the SCARGLE and CLEAN programs 
(Scargle 1982; Roberts et al. 1987) found a most likely periodicity in the range 2.5 to 3.7 
hours, which is close to the expected rotational period of LP 944— 20 (vsin i=30 kms -1 , 
P=4.4 hours sin i for R=Rj U pit er )- The CLEAN periodogram is shown in (Fig. 2). We put 
the UVES RV data in phase with a period of 3.7 hours (Fig. 3) and found a good fit with a 
sinusoid curve with an amplitude of 3.5 kms -1 . 

The NIRSPEC RV data, however, do not follow the periodic trend of the UVES data 
(Fig. 3), and hence the hypothesis of a planet to explain the UVES RV variations can be 
discarded because the RV signature of a planet should not depend on wavelength. In fact, our 
NIRSPEC data has an rms dispersion of 0.36 kms -1 , which is consistent with observational 
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uncertainties. Using a mass of 60 Jupiters for LP 944—20, we can rule out the presence of 
any planet with mass larger than m sini~1.0 Jovian masses and period less than 30 days 
using the NIRSPEC RV dataset. 

We confirm that the RV of LP 944—20 is variable at optical wavelengths as previously 
reported from UVES observations (Guenther & Wuchterl 2003). These authors also noted 
that the RV variability could not be due to starspots because the filling factor would have 
to be too large. Now that we have ruled out the existence of a giant planet, there seems 
to be only one possibility left, which is the presence of inhomogeneous surface features 
that are not like starspots. Theoretical models have predicted that there could be dust 
clouds in BDs, and several searches for photometric variability have been performed among 
late-M and L dwarfs (Bailer- Jones & Mundt 2001; Martin et al. 2001; Gelino et al. 2002; 
Goldman 2005). In fact, photometric variability has been detected in LP 944— 20 at opti- 
cal wavelengths (Tinney & Tolley 1999). On the other hand, near-infrared monitoring of 
about 20 late-M and L dwarfs failed to detect any variability (Bailer- Jones & Lamm 2003; 
Koen et al. 2005). 

LP944-20 and other ultracool dwarfs dramatically violate the X-ray/radio-emission 
relation for normal stars, as its quiescent radio emission is 4 to 5 orders of magnitude 
brighter than expected from its X-ray emission (Berger et al. 2005). The radio emission 
observed in some ultracool dwarfs imply a large magnetic field strength (over a hundred 
Gauss averaged over the whole surface). Thus, while it is now clear that at least some 
old BDs have strong magnetic fields, the properties of the resulting coronae are strik- 
ingly different from those of stars. Possibly, not only the coronae differ from those of 
stars but also the structure of the magnetic field itself. As shown by several authors 
(Chabrier & Kuker 2006; Dobler et al. 2006) fields of fully convective objects (like BDs) 
are expected not to be concentrated in small spots but to be distributed on a global scale. 
Additionally, due to the low temperature of old BDs the degree of ionisation in the at- 
mosphere is very low, leading to a very low degree of the coupling between the magnetic 
field and the atmosphere. The coupling between the gas and the magnetic field is usu- 
ally described in terms of the magnetic Reynolds number R m = lv/r] (were / is a length 
scale, v a velocity scale, and r\ the magnetic diffusivity (Priest 1982)). While for the sun 

R m is of the order of unity, for ultracool dwarfs R m could be as low as 10~ 15 10~ 20 

(Meyer & Meyer-Hofmeister 1999; Mohanty et al. 2002). If this is true, BDs with strong 
magnetic field should not have spots, and thus the gas-flow in the atmosphere should not be 
effected by the magnetic field. 

We propose here that the RV-variability detected with UVES may be due to weather 
effects, i.e. variability of cloud coverage and/or flows, streams or long-lived nonmagnetic 
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spots in the upper photosphere of this object. The near-infrared observations probe deeper in 
the photosphere, and hence they are less affected by those motions or structures. The small 
amplitude of the photometric and spectroscopic variability reported by Tinney & Trolley 
(1999) and Guenther & Wuchterl (2003), respectively, indicate that the surface features do 
not have a large contrast with the photosphere, and hence they are not like starspots. 

Independently from a deep physical understanding of the cause of the RV variations in 
the optical, we can safely conclude on pure empirical grounds that future RV-based planet 
searches in late-M stars and BDs will be better conducted at near-infrared wavelengths 
than in the optical because the targets are brighter and because the measurements are less 
affected by photospheric noise. This finding encourages the use of a new generation of high- 
resolution near-infrared spectrographs (Kaufl et al. 2004; Martin et al. 2005) for surveying 
the planetary population of very low-mass stars and BDs. 
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Fig. 1. — Average NIRSPEC spectra for our target (LP 944— 20) and the RV template 
(VB10). The spectra have been shifted to vacuum wavelength and are offset in the vertical 
axis. The upper panel shows a NIRSPEC order dominated by FeH molecular lines, whereas 
the lower panel shows an order which is dominated by a KI doublet and FeH molecular lines. 
Both of them were used to measure the radial velocity of this brown dwarf. 
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Fig. 2. — CLEAN periodogram analysis of our UVES RV dataset of LP 944— 20. 5 iterations 
were used and the gain was 0.1. The highest peak is marked. 
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Fig. 3. — RV curve of LP 944— 20 folded with a period of 3.7 hours found in the VLT 
data. A sinusoidal curve provides a good fit for the UVES data (filled circles) but not for 
the NIRSPEC data (open circles). The individual NIRSPEC points for each spectrum are 
shown, whereas nightly averages are given in Table 1. 
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Table 1. RV data for LP 944-20 



HJD-2400000.0 


RV {km/s) 


a (km/sj 


Instrument 


51890.80439 


9.04 


1.00 


NIRSPEC 


52189.74347 


9.09 


0.78 


UVES 


52213.62965 


9.18 


0.89 


UVES 


52214.60754 


11.42 


0.73 


UVES 


52224.77031 


8.96 


0.46 


UVES 


52236.73029 


12.86 


0.43 


UVES 


52237.83445 


8.33 


0.57 


UVES 


52915.87534 


10.58 


1.68 


UVES 


52920.63108 


10.51 


0.72 


UVES 


52927.80874 


11.61 


1.23 


UVES 


53006.74733 


7.35 


0.94 


UVES 


53010.58583 


7.89 


0.76 


UVES 


53010.70228 


5.05 


2.17 


UVES 


53017.64839 


9.00 


1.21 


UVES 


53022.58255 


9.37 


0.60 


UVES 


53030.63830 


8.16 


0.80 


UVES 


53669.93870 


8.25 


0.43 


NIRSPEC 


53670.97038 


8.69 


0.40 


NIRSPEC 


53671.96188 


8.61 


0.48 


NIRSPEC 


53672.94703 


9.50 


1.56 


NIRSPEC 


53753.70181 


8.70 


0.77 


NIRSPEC 



